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We argue that the electron star and the AdS Dirac hair solution are two limits of the free charged 
Fermi gas in AdS. Spectral functions of holographic duals to probe fermions in the background of 
electron stars have a free parameter that quantifies the number of constituent fermions that make 
up the charge and energy density characterizing the electron star solution. The strict electron star 
limit takes this number to be infinite. The Dirac hair solution is the limit where this number is 

(N 

£> unity. This is evident in the behavior of the distribution of holographically dual Fermi surfaces. As 

00 

we decrease the number of constituents in a fixed electron star background the number of Fermi 

re- 
surfaces also decreases. An improved holographic Fermi groundstate should be a configuration that 

shares the qualitative properties of both limits. 
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I. INTRODUCTION 



The insight provided by the application of the AdS / CFT correspondence to finite density 
Fermi systems has given brand new perspectives on the theoretical robustness of non- Fermi 
liquids [IH3]; on an understanding of the non-perturbative stability of the regular Fermi 
liquid equivalent to order parameter universality for bosons [U [5] , and most importantly on 
the notion of fermionic criticality: Fermi systems with no scale. In essence strongly coupled 
conformally invariant fermi systems are one answer to the grand theoretical question of 
fermionic condensed matter: Are there finite density Fermi systems that do not refer at any 
stage to an underlying perturbative Fermi gas? 

It is natural to ask to what extent AdS/CFT can provide a more complete answer to this 
question. Assuming, almost tautologically, that the underlying system is strongly coupled 
and there is in addition some notion of a large N limit, the Fermi system is dual to clas- 
sical general relativity with a negative cosmological constant coupled to charged fermions 
and electromagnetism. As AdS/CFT maps quantum numbers to quantum numbers, finite 
density configurations of the strongly coupled large N system correspond to solutions of this 
Einstein-Maxwell-Dirac theory with finite charge density. Since the AdS fermions are the 
only object carrying charge, and the gravity system is weakly coupled, one is immediately 
inclined to infer that the generic solution is a weakly coupled charged Fermi gas coupled to 
AdS gravity: in other words an AdS electron star [SJ E] , the charged equivalent of a neutron 
star in asymptotically anti-de Sitter space jS. 9j . 

Nothing can seem more straightforward. Given the total charge density Q of interest, 
one constructs the free fermionic wavefunctions in this system, and fills them one by one in 
increasing energy until the total charge equals Q. For macroscopic values of Q these fermions 
themselves will backreact on the geometry. One can compute this backreaction; it changes 
the potential for the free fermions at subleading order. Correcting the wavefunctions at 
this subleading order, one converges on the true solution order by order in the gravitational 
strength K 2 Ej ull system . Here Ef u u system is the energy carried by the Fermi system and k 2 is 
the gravitational coupling constant k 2 = 8tcG Newton in the AdS gravity system. Perturbation 
theory in k is dual to the 1/N expansion in the associated condensed matter system. 

The starting point of the backreaction computation is to follow Tolman-Oppenheimer- 
Volkov (TOV) and use a Thomas-Fermi (TF) approximation for the lowest order one-loop 



2 



contribution [6H9]. The Thomas- Fermi approximation applies when the number of con- 
stituent fermions making up the Fermi gas is infinite. For neutral fermions this equates to 
the statement that the energy-spacing between the levels is neglible compared to the chem- 
ical potential associated with Q, AE/fi — > 0. For charged fermions the Thomas-Fermi limit 
is more direct: it is the limit q/Q — > where q is the charge of each constituent fermion. 1 
This has been the guiding principle behind the approaches [6T4TT] and the recent papers 
[T2| [13] , with the natural assumption that all corrections beyond Thomas- Fermi are small 
quantitative changes rather than qualitative ones. On closer inspection, however, this com- 
pletely natural TF-electron star poses a number of puzzles. The most prominent perhaps 
arises from the AdS/CFT correspondence finding that every normalizable fermionic wave- 
function in the gravitational bulk corresponds to a fermionic quasiparticle excitation in the 
dual condensed matter system. In particular occupying a particular wavefunction is dual to 
having a particular Fermi-liquid state [I]. In the Thomas Fermi limit the gravity dual thus 
describes an infinity of Fermi liquids, whereas the generic condensed matter expectation 
would have been that a been that a single(/few) liquid(s) would be the generic groundstate 
away from the strongly coupled fermionic quantum critical point at zero charge density. This 
zoo of Fermi surfaces is already present in the grand canonical approaches at fixed fi (ex- 
tremal AdS-Reissner-Nordstrom (AdS-RN) black holes) [3] and a natural explanation would 
be that this is a large N effect. This idea, that the gravity theory is dual to a condensed 
matter system with N species of fermions, and increasing the charge density "populates" 
more and more of the distinct species of Fermi liquids, is very surprising from the condensed 
matter perspective. Away from criticality one would expect the generic groundstate to be a 
single Fermi-liquid or some broken state due to pairing. To pose the puzzle sharply, once one 
has a fermionic quasiparticle one should be able to adiabatically continue it to a free Fermi 
gas, which would imply that the free limit of the strongly coupled fermionic CFT is not a 
single but a system of order N fermions with an ordered distribution of fermi-momenta. A 
possible explanation of the multitude of Fermi surfaces that is consistent with a single Fermi 
surface at weak coupling is that AdS/CFT describes so-called "deconfined and/or fraction- 
alized Fermi-liquids" where the number of Fermi surfaces is directly tied to the coupling 

1 For a fermion in an harmonic oscillator potential E n — h(n— \/2)uj: thus AE/E tota i = l/J^f ( n ~ V^) = 
2 /N 2 . 
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strength [T2HT6] . It would argue that fermionic quantum criticality goes hand in hand with 
fractionalization for which there is currently scant experimental evidence. 

The second puzzle is more technical. Since quantum numbers in the gravity system equal 
the quantum numbers in the dual condensed matter system, one is inclined to infer that each 
subsequent AdS fermion wavefunction has incrementally higher energy than the previous one. 
Yet analyticity of the Dirac equation implies that all normalizable wavefunctions must have 
strictly vanishing energy [IT]. It poses the question how the order in which the fermions 
populate the Fermi gas is determined. 

The third puzzle is that in the Thomas Fermi limit the Fermi gas is gravitationally strictly 
confined to a bounded region: famously, the TOV-neutron star has an edge. In AdS/CFT, 
however, all information about the dual condensed matter system is read off at asymptotic 
AdS infinity. Qualitatively, one can think of AdS/CFT as an "experiment" analogous to 
probing a spatially confined Fermi gas with a tunneling microscope held to the exterior of 
the trap. Extracting the information of the dual condensed matter system is probing the 
AdS Dirac system confined by a gravitoelectric trap instead of a magneto-optical trap for 
cold atoms. Although the Thomas- Fermi limit should reliably capture the charge and energy 
densities in the system, its abrupt non-analytic change at the edge (in a trapped system) and 
effective absence of a density far away from the center are well known to cause qualitative 
deficiencies in the description of the system. Specifically Friedel oscillations — quantum 
interference in the outside tails of the charged fermion density, controlled by the ratio q/Q 
and measured by a tunneling microscope — are absent. Analogously, there could be qualita- 
tive features in the AdS asymptotics of both the gravito-electric background and the Dirac 
wavefunctions in that adjusted background that are missed by the TF-approximation. The 
AdS asymptotics in turn specify the physics of the dual condensed matter system and since 
our main interest is to use AdS / CFT to understand quantum critical fermion systems where 
q/Q is finite, the possibility of a qualitative change inherent in the Thomas Fermi limit 
should be considered. 

There is another candidate AdS description of the dual of a strongly coupled finite density 
Fermi system: the AdS black hole with Dirac hair [HE]. One arrives at this solution when one 
starts one's reasoning from the dual condensed matter system, rather than the Dirac fields 
in AdS gravity. Insisting that the system collapses to a generic single species Fermi-liquid 
ground state, the dual gravity description is that of an AdS Einstein-Dirac-Maxwell system 
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with a single nonzero normalizable Dirac wavefunction. To have a macroscopic backreaction 
the charge of this single Dirac field must be macroscopic. The intuitive way to view this 
solution is as the other simplest approximation to free Fermi gas coupled to gravity. What 
we mean is that the full gravito-electric response is in all cases controlled by the total charge 
Q of the solution: as charge is conserved it is proportional to the constituent charge q times 
the number of fermions np AdS and the two simple limits correspond to np — > oo,q — > 
with Q = qnp fixed or np — > 1, q — >■ Q. The former is the Thomas-Fermi electron star, 
the latter is the AdS Dirac hair solution. In the context of AdS / CFT there is a significant 
difference between the two solutions in that the Dirac Hair solution clearly does not give rise 
to the puzzles 1, 2 and 3: there is by construction no zoo of Fermi-surfaces and therefore no 
ordering. Moreover since the wavefunction is demanded to be normalizable, it manifestly 
encodes the properties of the system at the AdS boundary. On the other hand the AdS 
Dirac hair solution does pose the puzzle that under normal conditions the total charge Q is 
much larger than the constituent charge q both from the gravity /string theory point of view 
and the condensed matter perspective. Generically one would expect a Fermi gas electron 
star rather than Dirac hair. 

In this article we shall provide evidence for this point of view that the AdS electron star 
and the AdS Dirac hair solution are two limits of the same underlying system. Specifically 
we shall show that (1) the electron star solution indeed has the constituent charge as a 
free parameter which is formally sent to zero to obtain the Thomas-Fermi approximation. 
(2) The number of normalizable wavefunctions in the electron star depend on the value of 
the constituent charge q. We show this by computing the electron star spectral functions. 
They depend in similar way on q as the first AdS/CFT Fermi system studies in an AdS- 
RN background. In the formal limit where q — > Q, only one normalizable mode remains 
and the spectral function wavefunction resembles the Dirac Hair solution, underlining their 
underlying equivalence. Since both approximations have qualitative differences as a descrip- 
tion of the AdS dual to strongly coupled fermionic systems, it argues that an improved 
approximation which has characteristics of both is called for. 

The results here are complimentary to and share an analysis of electron star spectral 
functions with the two recent articles [12] and [13] that appeared in the course of this work 
(see also [IB] for fermion spectral functions in general Lifshitz backgrounds). Our motivation 
to probe the system away from the direct electron star limit differs: we have therefore been 
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more precise in defining this limit and in the analysis of the Dirac equation in the electron 
star background. 



II. EINSTEIN-MAXWELL THEORY COUPLED TO CHARGED FERMIONS 



The Lagrangian that describes both the electron star and Dirac Hair approximation is 
Einstein-Maxwell theory coupled to charged matter 

1 ,„ 6 \ f 



S 



2k 2 



(R+j^)-^- 2 F 2 + C matter (e^A,) 



(2-1) 



where L is the AdS radius, q is the electric charge and k is the gravitational coupling 
constant. It is useful to scale the electromagnetic interaction to be of the same order as the 
gravitational interaction and measure all lengths in terms of the AdS radius L: 



(2.2) 



The system then becomes 



/ 



2k 2 



R + 6-^F 2 ) + L 4 £ matter (L< q —A^) 



K 



(2.3) 



Note that in the rescaled variables the effective charge of charged matter now depends on 
the ratio of the electromagnetic to gravitational coupling constant: q eff = qL/n. For the 
case of interest, charged fermions, the Lagrangian in these variables is 



r A r (jp A ) - 

-"—f ermions V-^ ^ ? M/ 



Ap) - mL 



tt, (2.4) 



where ^ is defined as — i^/^T . Compared to the conventional normalization the Dirac 
field has been made dimensionless \I/ = ftV^V'conventionai- With this normalization all terms 
in the action have a factor L 2 / k 2 and it will therefore scale out of the equations of motion 



Rfjiu — ~j9iivR — 3^i/ — {^F iip F v p — -g^FpeF*" 7 + T f 
D^F^ = -q eii Jf ermions 



fermions \ 

pu 



J 5 



(2.5) 



with 



T, 



fermions 



9 ») + 



BC-p 
1 BC 



K 



K 2 L 2 



fermions i 



" fermions 



(2.6) 
(2.7) 
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where the symmetrization is defined as BiJJ v \ = B^C V + B U C^ and the Dirac equation 



e» A T A (d» + ]^ C T BC - i q —A,) - mL 



* = 0. 



(2.8) 



The stress-tensor and current are to be evaluated in the specific state of the system. 



For a single excited wavefunction, obeying (2.8), this gives the AdS Dirac hair solution 



constructed in [I]. (More specifically, the Dirac hair solution consists of a radially isotropic 
set of wavefunctions with identical momentum size \k\ = ^Jk\ + k^, such that the Pauli 
principle plays no role.) For multiple occupied fermion states, even without backreaction 



due to gravity, adding the contributions of each separate solution to (2.8) rapidly becomes 



very involved. In such a many-body-system, the collective effect of the multiple occupied 
fermion states is better captured in a "fluid" approximation 



rp fluid 



N, 



fluid 



nu. 



(2.9) 



with 



p (u^T^yU )matter only j H i^fiJ^ matter only 



(2.10) 



In the center-of-mass rest frame of the multiple fermion system (u^ = (eto, 0, 0, 0)), the 
expressions for the stress-tensor and charge density are given by the one-loop equal-time 
expectation values (as opposed to time-ordered correlation functions) 

1 



p = m^T^dt + -u? B ?ab - iq etl A t )^{t)). 



(2.11) 



By the optical theorem the expectation value is equal to twice imaginary part of the Feynman 



propagator 



p = lim 2ImTr 



elY^d t + -u? B T AB - iq eti A t )Gf b (t\ t) 



(2.12) 



In all situations of interest, all background fields will only have dependence on the radial 
AdS direction; in that case the spin connection can be absorbed in the normalization of the 
spinor wavefunction. 3 In an adiabatic approximation for the radial dependence of eto and A t 



2 From unitarity for the S matrix S — 1 one obtains the optical theorem T^T = 2ImT for the transition 
matrix T defined as 5* = 1 + iT. 

3 i.e. one can redefine spinors x(r) = /(r)^(r) such that the connection term is no longer present in the 
equation of motion. 
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- where yU l0C (r) = g eff eo(r)At(r) and cu(r) = — ieg(r)c^; — this yields the known expression 
for a many-body-fermion system at finite chemical potential 

d 3 kdu 



p(r) = lim 2 

f3— >oo 



lim 

/9— >oc 



Mr) - /Xi oc (r)] lmTrT^4(w, fc) 



(2vr) 4 
dkdu 



471-3 



[/c 2 (o; — /i 



1 1 T.// 3 / 

- - -tanh(-(o;-//)) 



Tr(*r°) 2 ^-7r5((u; - /i) - v 7 *: 2 + (mi) 2 ) 



lim — — 

/3^oo ir 2 L 2 

1 k 2 ^° 



duf FD ((3(u - fj)) [(u - ;u) 2 - (mi) 2 ] [w - /i] 
dEE 2 ^E 2 - (mi) 2 . 



(a> — fi)9(oj — jj, — mL) 
\f(w — fi) 2 — (mL) 2 



71-2 L 2 j mL 

The normalization k 2 / L 2 follows from the unconventional normalization of the Dirac field 



(2.13) 



in eq. (2.4). Similarly 



n 



^L 2 



Mloc 



dEE^E 2 -{mL) 2 = ±- 2 H -^ 

mL ^ 



loc 



(mL) 



2\3/2 



(2.14) 



The adiabatic approximation is valid for highly localized wavef unctions, i.e. the expres- 
sion must be dominated by high momenta (especially in the radial direction). The exact 
expression on the other hand will not have a continuum of solutions to the harmonic con- 
dition — T°uj + Pfcj + T z k z — T°fi\ oc — imL = 0. Normalizable solutions to the AdS Dirac 
equations only occur at discrete momenta — one can think of the gravitational background 
as a potential well. The adiabatic approximation is therefore equivalent to the Thomas- Fermi 
approximation for a Fermi-gas in a box. 

To get an estimate for the parameter range where the adiabatic approximation holds, 
consider the adiabatic bound d r fi loc (r) <C /Ui oc (r) 2 . Using the field equation for Aq = 

A t loc/9eff : 



d 2 [J>ioc 



7eff" ' 



(2.15) 



this bound is equivalent to requiring 



loc 



(—) 2 n < 2fi loc d r fi loc 

K 



(2.16) 



One can see this readily by converting the dimensionless definition of p, eq (2.11 ), to the standard dimen 



sion. Using capitals for dimensionless quantities and lower-case for dimensionful ones 



p ~ (*a T *) ~ k 2 L 2 (?MV>) ~ k 2 L 2 r dee 2 \fe 
with fiL — /ii oc above. 



V- 



^ dEE 2 \J E 2 - {mL) 2 



8 



where in the last line we used the original bound again. If the chemical potential scale 



is considerably higher than the mass of the fermion, we may use (2.14) to approximate 



n ~ f^ioc- Thus the adiabatic bound is equivalent to 

9eff K 



q = < 1 (2.17) 

the statement that the constituent charge of the fermions is infinitesimal. Note that in the 



^spatial AdS ^Kspatial AdS ^^spatial AdS 

1 f , i 9eff« 



rescaled action (2.3, 2.4), L/k plays the role of 1/H, and eq. (2.17) is thus equivalent to the 
semiclassical limit h—>0 with q eii fixed. Since AdS/CFT relates L/k ~ N c this acquires the 
meaning in the context of holography that there is a large N c scaling limit [T2l ITB"] of the CFT 
with fermionic operators where the RG-flow is "adiabatic" . Returning to the gravitational 
description the additional assumption that the chemical potential is much larger than the 
mass is equivalent to 

/Ototal r/Qtotal r r 

" phy8 dr^^i^n) (2.18) 



^spatial AdS 

This implies that the total charge density in AdS is much larger than that of a single charged 
particle (as long as mL ~ 1). The adiabatic limit is therefore equivalent to a thermodynamic 
limit where the Fermi gas consists of an infinite number of constituents, n — > oo, q — > such 
that the total charge Q ~ nq remains finite. 

The adiabatic limit of a many-body fermion system coupled to gravity are the Tolman- 
Oppenheimer-Volkov equations. Solving this in asymptotically AdS gives us the charged 
neutron or electron star constructed in . Knowing the quantitative form of the adiabatic 
limit, it is now easy to distinguish the electron star solution from the "single wavefunction" 
Dirac Hair solution. The latter is trivially the single particle limit n — > 1, q — y Q with 
the total charge Q finite. The electron star and Dirac Hair black hole are opposing limit- 
solutions of the same system. We shall now make this connection more visible by identifying 
a formal dialing parameter that interpolates between the two solutions. 

To do so we shall need the full adiabatic Tolman-Oppenheimer-Volkov equations for the 
AdS electron star [7]. Since the fluid is homogeneous and isotropic, the background metric 
and electrostatic potential will respect these symmetries and will be of the form (recall that 



we are already using "dimensionless" lengths, eq. (2.2)) 

ds 2 = -f(r)dt 2 + g(r)dr 2 + r 2 (dx 2 + dy 2 ), A = h(r)dt, (2.19) 



where f{r),g(r), h(r) are functions of r; the horizon is located at r = and the boundary 
is at r = oo. Combining this ansatz with a rescaling mL = q eii rn the bosonic background 
equations of motion become [7] 



1 f f' 9'\ gha q\ tt K 2 f VJ 



+ 777 ~~ 9(3 + P) + -7, = 0, cr = / deeVe 2 - m 2 , 



rf 2f * K " r* ' tt 2 L 2 
, „ 2 , . qo~ ( rhh' „\ /i , 

ft + r ft -^(^ + / j= ' " P = "-77 ff ' (2 ' 20) 

where we have used that fi loc = q e tth/y/J and a = nq efi is the rescaled local charge 
density. What one immediately notes is that the Tolman-Oppenheimer-Volkov equations of 
motion for the background only depend on the parameters p = *" L 2 and m, whereas the 
original Lagrangian and the fermion equation of motion also depend on q eft - ' 



It is therefore natural to guess that the parameter q Btt = qL/n will be the interpolating 
parameter away from the adiabatic electron star limit towards the Dirac Hair BH. 



Indeed in these natural electron star variables the adiabatic bound (2.17) translates into 



$ < ^ = % ■ (2.21) 



Thus we see that for a given electron star background with (3 fixed decreasing k/L improves 
the adiabatic fluid approximation whereas increasing k/L makes the adiabatic approximation 
poorer and poorer. "Dialing k/L up/down" therefore interpolates between the electron 
star and the Dirac Hair BH. Counterintuively improving adiabaticity by decreasing k/L 
corresponds to increasing q etf for fixed q, but this is just a consequence of recasting the 
system in natural electron star variables. A better way to view improving adiabaticity is to 
decrease the microscopic charge q but while keeping q ett fixed; this shows that a better way 
to think of q ei ± is as the total charge rather than the efffective constituent charge. 

The parameter k/L = q/q e ±± parametrizes the gravitational coupling strength in units 
of the AdS curvature, and one might worry that "dialing k/L up" pushes one outside the 
regime of classical gravity. This is not the case. One can easily have $ 3> 1 and tune k/L 
towards or away from the adiabatic limit within the regime of classical gravity. From eq. 



(2.17) we see that the edge of validity of the adiabatic regime $ ~ L 2 /k 2 is simply equivalent 
to a microscopic charge q = 1 which clearly has a classical gravity description. It is not hard 
to see that the statement above is the equivalent of changing the level splitting in the Fermi 
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gas, while keeping the overall energy/charge fixed. In a Fermi gas microscopically both the 
overall energy and the level splitting depends on h. Naively increasing h increases both, but 
one can move away from the adiabatic limit either by decreasing the overall charge density, 
keeping h fixed or by keeping the charge density fixed and raising h. Using again the analogy 
between k/L and h, the electron star situation is qualitatively the same where one should 
think of $ ~ q A L 2 / k 2 parametrizing the microscopic charge. One can either insist on keeping 
k/L fixed and increase the microscopic charge (3 to increase the level splitting or one can 
keep (3 fixed and increase k/L. In the electron star, however, the background geometry 
changes with (3 in addition to the level splitting, and it is therefore more straightforward to 
keep /3 and the geometry fixed, while dialing k/L. 

We will now give evidence for our claim that the electron star and Dirac Hair solution are 
two opposing limits. To do so, we need to identify an observable that goes either beyond the 
adiabatic background approximation or beyond the single particle approximation. Since the 
generic intermediate state is still a many-body fermion system, the more natural starting 
point is the electron star background and perturb away from there. Realizing then that 
the fermion equation of motion already depends directly on the dialing parameter q eti the 
obvious observables are the single fermion spectral functions in the electron star background. 
Since one must specify a value for q eii to compute these, they directly probe the microscopic 
charge of the fermion and are thus always beyond the strict electron star limit q — > 0. In the 
next two sections we will compute these and show that they indeed reflect the interpretation 
of <7ef f as the interpolating parameter between the electron star and Dirac Hair BH. 



III. FERMION SPECTRAL FUNCTIONS IN THE ELECTRON STAR BACK- 
GROUND 



To compute the fermion spectral functions in the electron star background we shall choose 
a specific representative of the family of electron stars parametrized by $ and m. Rather 
than using (3 and 771 the metric of an electron star is more conveniently characterized by its 



Lifshitz-scaling behavior near the interior horizon r — > 0. From the field equations (2.20) 
the limiting interior behavior of f(r),g(r), h(r) is 

f(r) = r 2z + . . . , g(r) = ^ + . . . , h(r) = h^r* + ... (3.1) 
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The scaling behavior is determined by the dynamical critical exponent z, which is a function 
of (3, rh [7] and it is conventionally used to classify the metric instead of (3. The full 
electron star metric is then generated from this horizon scaling behavior by integrating up 
an irrelevant RG-flow [HI [20] 

/ = r 2 ^l + Ar- Q + ...^, g = ^(l + g ir -* + ...^ h = h 00 r z (l + h 1 r- a + ...y.2) 
with 

2 + z J9z s - 21z 2 + 40z - 28 - m 2 zU - 3z) 2 . . 

a = . . (3.3) 

2 2y/(l -m 2 )z- 1 

Scaling fx — > bfi is equal to a coordinate transformation r — > b l ' a r and t — > b z ^ a t, and the 
sign of fi is fixed to be negative in order to be able to match onto an asymptotically AdS4 
solution. Thus fi = —1 and gi and hi are then uniquely determined by the equations of 
motion. 

Famously, integrating the equations of motion up the RG-flow outwards towards the 
boundary fails at a finite distance r s . This is the edge of the electron star. Beyond the edge 
of the electron star, there is no fluid present and the spacetime is that of an AdS4-RN black 
hole with the metric 

, 22 M Q 2 c 2 , Q 

f = cr + #=7' h = ^--- ( 3 - 4 ) 

Demanding the full metric is smooth at the radius of electron star r s determines the constants 
c, M and Q. The dual field theory is defined on the plane ds 2 = —c 2 dt 2 + dx 2 + dy 2 . 

The specific electron star background we shall choose without loss of generality is the 
one with z = 2,m = 0.36 (Fig. [T]) 5 , smoothly matched at r s ~ 4.25252 onto a AdS-RN 
black- hole. 

The CFT fermion spectral functions now follow from solving the Dirac equation in this 
background [U |2] 



e^r A ( df, + ^-c^ABr^ 8 - iq etf ) - m eff 



^ = (3.5) 



5 This background has c ~ 1.021, M ~ 3.601, Q ~ 2.534, /i ~ 2.132, $ ~ 19.951, ^ ~ 1.887, ft 
l/v^, a ^ -1.626, /i = -1, g x ~ -0.4457, fti ~ -0.6445. 
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FIG. 1: Electron star metric for z = 2,m = 0.36, c ~ 1.021, M ~ 3.601, Q ~ 2.534, (t ~ 2.132 
compared to pure AdS. Shown are f(r)/r 2 (Blue), r 2 g(r) (Red) and h(r) (Orange). The asymptotic 
AdS-RN value of h(r) is the dashed blue line. For future use we have also given fi loc = h/yff 
(Green) and fj, qeft = y/g™h/y/f (Red Dashed) At the edge of the star r s ~ 4.253 (the intersection 
of the purple dashed line setting the value of m eff with fii oc ) one sees the convergence to pure AdS 
in the constant asymptotes of f(r)/r 2 and r 2 g(r). 

where q eii and m eif in terms of the parameters of the electron star equal 

?eff = I k2 I , meff = q eff m = m I ^ I . (3.6) 

For a given electron star background, i.e. a fixed /3, rh the fermion spectral function will 
therefore depend on the ratio L/k. For L/k ^> (3 1 ^ 2 the poles in these spectral functions 
characterize the occupied states in a many-body gravitational Fermi system that is well 
approximated by the electron star. As L/k is lowered for fixed j3 the electron star background 
becomes a poorer and poorer approximation to the true state and we should see this reflected 
in both the number of poles in the spectral function and their location. 
Projecting the Dirac equation onto two-component T- eigenspinors 

and using isotropy to set k y = 0, one can choose a basis of Dirac matrices where one obtains 
two decoupled sets of two simple coupled equations pQ 

V9u9 rr (d r =F m elf ^/g^)y ± = Ti(k x ~u)z T , (3.8) 

±i{k x + u)y± (3.9) 
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where u 



[uj + q c ffh). In this basis of Dirac matrices the CFT Green's function 



G = (Of + ij O^ + ) equals 



G = lim e 

£->0 



-2mL 




where £_) 



W- 



IZ- 

y+ 



(3.10) 



Rather than solving the coupled equations (3.8) it is convenient to solve for £± directly [TJ, 




d r C± = -2m eff v /^"^ ± =f (k x Tu)± (k x ± u)£±. 



(3.11) 

For the spectral function A = ImTrG^ we are interested in the retarded Green function. 
This is obtained by imposing in-falling boundary conditions near the horizon r = 0. Since 
the electron star is a "zero-temperature" solution this requires a more careful analysis than 
for a generic horizon. To ensure that the numerical integration we shall perform to obtain 



the full spectral function has the right infalling boundary conditions, we first solve eq. (3.11) 



to first subleading order around r = 0. There are two distinct branches. When u^O and 
k x r / u , r 2 / uj is small, the in-falling boundary condition near the horizon r = is (for z — 2) 



. k x r (k 2 x - 2im eit uj)r 
uj 2uj 2 

f_( r ) = i + i + — J — + i 

uj 2uj 1 



l h 

2uj 
hk x r l ~< 



When uj — 0, i.e. k x r/uj is large, and r/k x — > 



_ 1 (geff^oo + m eff )r 



1 + 



kj. 

- m ef i)r 



UJ 

k-rT 



2uj 



UJ 



+ . 



(3.12) 



UJ 



2y/ 9ook x 
UJ 

2\/9oak x 



+ • 



+ 



(3.13) 



the boundary conditions (3.13) become real. As (3.11) are real equations, the spectral 



function vanishes in this case. This is essentially the statement that all poles in the Green's 
function occur at uj = [T7j. Note that the fact that the electron star uj = boundary 



conditions (3.11 ) are real for all values of k is qualitatively different from the AdS-RN uj = 



boundary conditions (eq. (26) in pQ). In the AdS-RN "quantum-critical" infrared governed 
by the near horizon AdS2 x M2 geometry, in general there is a special scale k Q below which 
the boundary condition turns complex. This scale k a is related to the surprising existence of 
an oscillatory region in the spectral function. One therefore infers that in a scale-ful Lifshitz 
infrared this oscillatory region is no longer present [131 EI]- We will confirm this in section 

El 
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FIG. 2: Electron star MDF spectral functions with multiple peaks as a function of k for u 



10 



-5 



2, rh = 0.36. The blue curve is for k = 0.091; the red curve is for k = 0.090. Note that 



the vertical axis is logarithmic. Visible is the rapidly decreasing spectral weight and increasingly 
narrower width for each successive peak as kp increases. 

A. Numerical results and discussion 

We can now solve for the spectral functions numerically. In Fig. [3] we plot the momentum- 
distribution-function (MDF) (the spectral function as a function of k) for fixed u = 10~ 5 , z = 
2, rh — 0.36 while changing the value of k. Before we comment on the dependence on 
q eff ~ kT 1 ! 2 which studies the deviation away from the adiabatic limit of a given electron 
star background (i.e. fixed dimensionless charge and fixed dimensionless energy density), 
there are several striking features that are immediately apparent: 

• As expected, there is a multitude of Fermi surfaces. They have very narrow width 
and their spectral weight decreases rapidly for each higher Fermi- momentum kp (Fig. 



2). This agrees with the exponential width T ~ exp(— (^) 



k ^l/(z-l), 



predicted by 



for gravitational backgrounds that are Lifshitz in the deep interior, which is the case 
for the electron star. This prediction is confirmed in [121 HH1 HE] and the latter two 
articles also show that the weight decreases in a corresponding exponential fashion. 
This exponential reduction of both the width and the weight as kp increases explains 
why we only see a finite number of peaks, though we expect a very large number. In 
the next section we will be able to count the number of peaks, even though we cannot 
resolve them all numerically. 
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• The generic value of kp of the peaks with visible spectral weight is much smaller than 
the effective chemical potential \i in the boundary field theory. This is quite different 
from the RN-AdS case where the Fermi momentum and chemical potential are of the 
same order. A numerical study cannot answer this, but the recent article p3] explains 
this. 6 

• Consistent with the boundary value analysis, there is no evidence of an oscillatory 
region. 

The most relevant property of the spectral functions for our question is that as k is increased 
the peak location kp decreases orderly and peaks disappear at various threshold values of k. 
This is the support for our argument that changing k changes the number of microscopic 
constituents in the electron star. Comparing the the behavior of the various Fermi momenta 
/cp in the electron star with the results in the extremal AdS-RN black-hole, they are quali- 
tatively identical when one equates n~ l l 2 ~ q eff with the charge of the probe fermion. We 
may therefore infer from our detailed understanding of the behavior of kp for AdS-RN that 
also for the electron star as kp is lowered peaks truly disappear from the spectrum until by 
extrapolation ultimately one remains: this is the AdS Dirac hair solution [3]. 

We can only make this inference qualitatively as the rapid decrease in spectral weight of 
each successive peak prevents an exact counting of Fermi surfaces in the numerical results 
for the electron star spectral functions. One aspect that we can already see is that as k 
decreases all present peaks shift to higher k, while new peaks emerging from the left for 
smaller kappa. This suggests a fermionic version of the UV/IR correspondence where the 
peak with lowest kp corresponds to the last occupied level, i.e. highest "energy" in the AdS 
electron star. We will now address both of these points in more detail. 

IV. FERMI SURFACE ORDERING: k F FROM A SCHRODINGER FORMULA- 
TION 

Our analysis of the behavior of boundary spectral functions as a function of k relies on the 
numerically quite evident peaks. Stricly speaking, however, we have not shown that there is 

6 In view of the verification of the Luttinger count for electron star spectra in [T2J [T3] , this had to be so. 
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(A) 




(Bl) k*io 3 (B2) k 

FIG. 3: (A) Electron-star MDF spectral functions as a function of k for z = 2,m = 0.36, u = 
10~ 5 . Because the peak height and weights decrease exponentially, we present the adjacent ranges 
k G [0.017,0.019] and k G [0.019,0.021] in two different plots with different vertical scale. (B1/B2) 
Locations of peaks of spectral functions as a function of k: comparison between the electron star 
(Bl) for z = 2, m = 0.36, u = 10~ 5 (the dashed gray line denotes the artificial separation in the 
3D representations in (A)) and AdS-RN (B2) for m = as a function of q in units where \i = \/3 
These two Fermi-surface 'spectra' are qualitatively similar. 

a true singularity in the Green's function at u; = 0,k = kp. We will do so by showing that the 
AdS Dirac equation, when recast as a Schrodinger problem has quasi-normalizable solutions 
at uj = for various k. As is well known, in AdS/CFT each such solution corresponds 
to a true pole in the boundary Green's function. Using a WKB approximation for this 
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Schrodinger problem we will in addition be able to estimate the number of poles for a fixed 
k and thereby provide a quantitative value for the deviation from the adiabatic background. 

We wish to emphasize that the analysis here is general and captures the behavior of 
spectral functions in all spherically symmetric and static backgrounds backgrounds alike, 
whether AdS-RN, Dirac hair or electron star. 



The cu = Dirac equation (3.5) for one set of components (3.8 3.9) with the replacement 
iy_ —7- equals 



V9u9 rr d r y- + meity/guV- = ~(k ~ frq eli )z+, 



V9ii9 rr d r z + - m eff y/g~iZ + = -(k + /t 9eff )?/_, (4.1) 

where p, qeli = J ^q e ffA t and we will drop the subscript x on k x . In our conventions z + 
(and y + ) is the fundamental component dual to the source of the fermionic operator in 
the CFT [TJ [2] . Rewriting the coupled first order Dirac equations as a single second order 
equation for z + \ 

d 2 r z + + Vd r z + + Qz + = 0, 
d r {9u9 rr ) d r fi q 



V 



2g ll g rr k + fi q> 



eff _ 



'eff 



rn eti d r Jg~ii rn e f fy /g^d r fi qaff 2 k 2 fi 2 

Q = 7 = + r~ —~ rn eff g rr —. 4.2 

y/9u9 rr k + n Qatt g u g rr 



the first thing one notes is that both V and Q diverge at some r = r* where fi qeff + k = 0. 
Since fi qeft is (chosen to be) a positive semidefinite function which increases from (x qeif = at 
the horizon, this implies that for negative k (with —k< fi qeff \ OQ ) the wavefunction is qualita- 
tively different from the wavefunction with positive k which experiences no singularity. The 
analysis is straightforward if we transform the first derivative away and recast it in the form 
of a Schrodinger equation by redefining the radial coordinate: 

- — °/>w (4 - 3) 

where Cq is an integration constant whose natural scale is of order cq ~ q~ fi . This is a 
simpler version of the generalized fc-dependent tortoise coordinate introduced in [3]. In the 



new coordinates the equation (4.2) is of the standard form: 

d 2 z + - V{s)z + = (4.4) 
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with potential 



-,fr 



<?ef f I 



The above potential (4.5) can also be written as 
1 



V(s) 



4( k + ^eJ 2 



{k 2 + ml if gu - ft 2 ) + m eff gii^fcrd r In v '"' 



(4.6) 



We note again the potential singularity for negative k, but before we discuss this we first 
need the boundary conditions. The universal boundary behavior is at spatial infinity and 
follows from the asymptotic AdS geometry. In the adapted coordinates r — > oo corresponds 
to s as follows from ds/dr ~ c (k + /t 9eff |oo)/ r2 - The potential therefore equals 

V{s) ~ — (m eff + m y + ... (4.7) 

and the asymptotic behavior of the two independent solutions equals z + = ais~ mett + 
^ iS i+m, ef f _)_____ second solution is normalizable and we thus demand a\ — 0. 
In the interior, the near-horizon geometry generically is Lifshitz 

ds 2 = -r 2z dt 2 + \dr 2 + r 2 (dx 2 + dy 2 ) + . . . , A = h^dt + . . . , (4.8) 

with finite dynamical critical exponent z — AdS-RN, which can be viewed as a special 
case where z — > oo, will be given separately. In adapted coordinates the interior r — > 
corresponds to s — > — oo and it is easy to show that in this limit potential behaves as 

V{s) ~ \ + \ (m eify /g^ + mefffl-oc - h^q^g^) + .... (4.9) 

Near the horizon the two independent solutions for the wavefunction z + therefore behave as 

z + -± a e- s/co + b e s/co . (4.10) 

The decaying solution a = is the normalizable solution we seek. 

Let us now address the possible singular behavior for k < 0. To understand what happens, 
let us first analyze the potential qualitatively for positive k. Since the potential is positive 



semi-definite at the horizon and the boundary, the Schrodinger system (4.4) only has a 



zero-energy normalizable solution if V(s) has a range s\ < s < S2 where it is negative. 



This can only at locations where k 2 < /t^ ff — m 2 it ga — m^gu^/g^dr In fc _^" ■ Defining 
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a "renormalized" position dependent mass m 2 = m 2 ff gn + m eff gaJg^d r In , this is 
the intuitive statement that the momenta must be smaller than the local chemical potential 
k 2 < fi 2 — m 2 eu . For positive k the saturation of this bound k 2 = ft 2 eff — m 2 en has at most 
two solutions, which are regular zeroes of the potential. This follows from the fact that 
fi 2 efi decreases from the boundary towards the interior. If the magnitude \k\ is too large 
the inequality cannot be satisfied, the potential is strictly positive and no solution exists. 
For negative k, however, the potential has in addition a triple pole at k 2 = fi 2 efl ; two poles 
arise from the prefactor and the third from the m ei fd r hi(k + /% ff ) term. This pole always 
occurs closer to the horizon than the zeroes and the potential therefore qualitatively looks 
like that in Fig. [1] (Since fi qetf decreases as we move inward from the boundary, starting with 
fi 2 eff > (x 2 qeii — yU 2 > k 2 , one first saturates the inequality that gives the zero in the potential 
as one moves inward.) Such a potential cannot support a zero-energy bound state, i.e. eq. 



(4.4) has no solution for negative k. In the case m efi = a double zero changes the triple 
pole to a single pole and the argument still holds. This does not mean that there are no 
k < poles in the CFT spectral function. They arise from the other physical polarization 
y + of the bulk fermion From the second set of decoupled first order equations for the 
other components of the Dirac equation (after replacing iz_ — > z_,) 

y/gug rr d r y + - m eii yfg~iy + — —(k — ft^)z-, 

V9ii9 rr d r z- + m eff y/g~iZ_ = —(k + fiq ell )y+, (4.11) 

and the associated second order differential EOM for y + : 

d 2 y + + Vd r y + + Q = 0, 
v _ d r {gng rr ) drAfcM 



2gug rr -k + fi q , 



eff 



rn &ii d r Jg~ii rn ef fy/g^d r fi qell 2 k p, 

Q = / = + 7- —~ m eff g rr 4.12 

V9ii9 rr -k + v Qeff gug rr 



one sees that the Schrodinger equation for y + is the k — > —k image of the equation (4.4) 
for z + and thus y + will only have zero-energy solutions for k < 0. For simplicity we will 
only analyze the z + case. Note that this semi-positive definite momentum structure of the 
poles is a feature of any AdS-to-Lifshitz metric different from AdS-RN, where one can have 
negative k solutions [3]. 



The exact solution of (4.4) with the above boundary conditions corresponding to poles 



in the CFT spectral function is difficult to find. By construction the system is however 
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V(s) 




FIG. 4: The behavior of the Schrodinger potential V(s) for z + when k is negative. Such a potential 
has no zero-energy bound state. The potential is rescaled to fit on a finite range. As \k\ is lowered 
below k max for which the potential is strictly positive, a triple pole appears which moves towards 
the horizon on the left (Fig A. The Blue, Red, Orange, Green curves are decreasing in |A;|). The pole 
hits the horizon for k = and disappears. Fig B. shows the special case m e ± f = where two zeroes 
collide with two of the triple poles to form a single pole. 



equivalent to a Schrodinger problem of finding a zero energy solution z + in the potential (4.5 ) 



and can be solved in the WKB approximation (see e.g. [31 123])- The WKB approximation 
holds when \d s V\ <C l^l 3 / 2 . Notice that this is more general than the background adiabacity 
limit m eff ^> l,q e t± 1 with /3, rh fixed . Combining background adiabaticity with a scaling 
limit k ^> 1, m e ±f ^> 1, q eif ^> 1 with c$k fixed and k is parametricaliy larger than fi qetl one 
recovers the WKB potential solved in [TJ1 US]- As our aim is to study the the deviation 
away from the background adiabatic limit we will be more general and study the WKB 
limit of the potential itself, without direct constraints on q eif , m ei± . And rather than testing 
the inequality \d s V\ <C | | 3 / 2 directly, we will rely on the rule of thumb that the WKB 
limit is justified when the number of nodes in the wave-function is large. We will therefore 
estimate the number n of bound states and use n > 1 as an empirical justification of our 
approach. 7 With this criterion we will be able to study the normalizable solutions to the 
Dirac equation/pole structure of the CFT spectral functions as a function of k/L. 

The potential is bounded both in the AdS boundary and at the horizon, and decreases 
towards intermediate values of r. We therefore have a standard WKB solution consisting of 



7 A large number of bound states n implies |<9 S V^| "C |V^| 3 / 2 if the potential has a single minimum, but as 
is well known there are systems, e.g. the harmonic oscillator, where the WKB approximation holds for 
small n as well. 
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three regions: 

• In the regions where V > 0, the solution decays exponentially: 

z + = c li2 \/- 1/4 exp(±^ dr'[c ^/g^ r (k + ft gatf )Vv]y (4.13) 

Here r±, r 2 are the turning points where V(ri) = = V(r 2 ). 

• In the region T\ < r < r 2 , i.e. V < 0, the solution is 



Z+ = C3 (-y)-V4 Re 



cxp 



(4.14) 



with the constant phase —in/ 4 originating in the connection formula at the turning 
point r 1 . 

Finding a WKB solution shows us that the peaks seen numerically are true poles in the 
spectral function. But it also allows us to estimate the number of peaks that the numerical 
approach could not resolve. The WKB quantization condition 



m r ' 

/ dr' c v 7 9 ii 9rr ( k + p, qetl ) V-V 

■J ri L 



7r(n + l/2) (4.15) 



counts the number of bound states with negative semi-definite energy. Note that n does not 
depend on the integral constant as there is also a factor l/c in y/—V. Since V depends on 
k, we will see that as we increase k this number decreases. The natural interpretation in 
the context of a bulk many-body Fermi system is that this establishes the ordering of the 
the filling of all the to = momentum shells in the electron star. For a fixed k one counts 
the modes that have been previously occupied and, consistent with our earlier deduction, 
the lowest /highest kp corresponds to the last/first occupied state. Though counterintuitive 
from a field theory perspective where normally E ~ kp , this UV/IR correspondence is very 
natural from the AdS-bulk, if one thinks of the electron star as a trapped electron gas. The 
last occupied state should then be the outermost state from the center, but this state has 
the lowest effective chemical potential and hence lowest kp. 

Let us now show this explicitly by analyzing the potential and the bound states in the 
electron star and AdS-RN. 
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(A) s ( B ) 

FIG. 5: The Schrodinger potential V(s) for the fermion component z + of in the ES background 
m = 0.36, z = 2, Co = 0.1. Fig. A. shows the dependence on the momentum k = 0.0185 (Purple), 
k = 5 (Blue), k = 10 (Red) for k = 0.092. Fig. B. shows the dependence on k = 0.086 (Purple), 
« = 0.092 (Blue), k = 0.1 (Red) for k = 0.0185. Recall that s = is the AdS boundary and 
s = — oo is the near-horizon region. 



Electron star 



The potential (4.6) for the electron star is given in Fig. |5jand the number of bound states 
as a function of k in Fig. |6j As stated the number of states decreases with increasing k, 
consistent with the analogy of the pole distribution of the spectral functions compared with 
AdS-RN. Moreover, we clearly see the significant increase in the number of states as we 
decrease k/L thereby improving the adiabaticity of the background. This vividly illustrates 
that the adiabatic limit corresponds to a large number of constituents. As all numbers of 
states are far larger than one, the use of the WKB is justified. 



The Reissner-Nordstrom case 

For AdS-RN the Schrodinger analysis requires a separate discussion of the near horizon 
boundary conditions, which we present here for completeness and comparison. Part of this 
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FIG. 6: The WKB estimate of the number of bound states n as a function of the momentum k for 
k = 0.086(Purple), 0.092(Blue), O.l(Red) (Fig A.); for k = 0.001 (Purple), 0.002(Blue), 0.003(Red) 
(Fig B.) and for k = 10~ 5 (Purple), 3 x 10" 5 (Blue), 5 x 10~ 5 (Red) (Fig C). Note the parametric 
increase in number of states as the adiabaticity of the background improves for smaller k. Both 
figures are for the electron star background with m = 0.36, z = 2. Since n 3> 1 in all cases, WKB 
gives a valid estimate. 



analysis is originally worked out in [3]. The AdS-RN black hole with metric 

ds 2 = L 2 (-f(r)dt 2 + ^ + r 2 (dx 2 + dy 2 )^ , (4.16) 

f(r) =r 2 (l + ^-^), (4.17) 

A = J I - -)dt, (4.18) 



r 

has near horizon geometry AdS2 x If 



dv 2 

ds 2 = -6(r - ifdt 2 + -, + (dx 2 + dy 2 ), (4.19) 

b(r — l) z 

A = y/str-Adt. (4.20) 



A coordinate redefinition of r in eq. (4.8) to r = (r^dSt ~ l) 1//2: shows that this corresponds 



to a dynamical critical exponent z = oo and is outside the validity of the previous analysis. 

Before we proceed, recall that the existence of AdS2 x IR 2 near-horizon region allows for a 
semi-analytic determination of the fermion spectral functions with the self-energy S ~ u 2uk F 
controlled by the IR conformal dimension 5k = 1/2 + i/& with 



When ujt is imaginary, which for q 2 > 2m 2 always happens for small k, the spectral function 
exhibits oscillatory behavior, but generically has finite weight at oj — 0. When is real, 
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FIG. 7: The Schrodinger potential V(s) for the fermion component z + of in the AdS-RN background 
r + = = \/3, gp = l,mL = 0.4, cq = 0.1. Fig. A. shows the dependence on the momentum 
k = 1 (Red), k = 2 (Purple), k = 3 (Blue) for charge g = 2.5. Fig. B. shows the dependence on the 
charge q — analogous to k in the ES background — . Shown are the values q = 2 (Blue), q = 2.5 
(Purple), q = 3 (Red) for the momentum k = 2. In both figures the Red potentials correspond to 
the oscillatory region u? < 0, the Purple potentials show the generic shape that can support an 
to = bound state, and the Blue potentials are strictly positive and no zero-energy bound state is 
present. Recall that s = is the AdS boundary and s = — oo is the near-horizon region. 



there are poles in the spectral functions at a finite number of different Fermi momenta 
kp. The associated quasiparticles can characterize a non-FL {y^ F < 1/2), a marginal FL 
(uk F = 1/2) or irregular FL (i>k F > 1/2) with linear dispersion but width r^w 2 [3]. 

The analytic form of the near-horizon metric allows us to solve exactly for the near 
horizon potential V in terms of s = ^%{k + q/y/2)\n(r — 1) + . . .. As noted in [3 J one 
remarkably obtains that the near-horizon potential for s — > — oo is proportional to the 
self-energy exponent: 

V{*)* 2{ ,J / ^ l + ---- (4-22) 
c 2 (k + q/V2) 2 

One immediately recognizes the oscillatory region v\ < of the spectral function as an 
to = Schrodinger potential which is "free" at the horizon s = — oo (Fig. [7]) and no bound 
state can form. Comparing with our previous results, we see that this oscillatory region is 
a distinct property of AdS-RN. For any Lifshitz near-horizon metric the potential is always 
positive-definite near the horizon and all to = solutions will be bounded, (see also [T2l IT3] ). 
As we increase k, v\ becomes positive, then the AdS-RN potential is also positive at the 
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FIG. 8: The qualitative behavior for negative k of the Schrodinger potential V(s) for the fermion 
component z + of in the AdS-RN background r + = = \/3,gp = l,mL = 0.1. The radial 
coordinate has been rescaled to a finite domain such that the full potential can be represented in 
the figure; on the right is the AdS boundary and left is the near-horizon region and the range is 
slightly extended beyond the true horizon, which is exactly at the short vertical line-segments on 
the right. Potentials are given for q = 12/\/3, k = —15 (Blue) for which the potential is strictly 
positive, k = — 10 (Red), k = —7 (Orange), which both have triple poles and the pole can be seen 
to move towards the horizon on the left as k decreases, and k = —4 (Green) which has no pole 
and a finite negative value at the horizon. The pole disappears for \k\ < q/V2 leaving a regular 
bounded potential which can support zero-energy bound states. 

horizon and bound zero-energy states can form. Increasing k further, one reaches a maximal 
kmax, above which the potential is always positive and no zero-energy bound state exists 
anymore. 

Because the near-horizon boundary conditions for AdS-RN differ from the general analy- 
sis, the possible singularity in the potential for k < also requires a separate study. This is 
clearly intimately tied to the existence of an oscillatory regime in the spectral function, as 
the previous analysis does apply for v\ > 0. The clearest way to understand what happens 
for v\ < is to analyze the potential explicitly. Again if \k\ > k max the potential is strictly 
positive definite, and no zero-energy bound state exists. As we decrease the magnitude of 
k < 0, a triple pole will form near the boundary when k = —fi qail (s), soon followed by a zero 
at k — — a/ £hq eU (s) 2 — m ren (s) 2 (see Fig. [4]). As we approach the horizon, in the general case 
where lim r _;.o /Aj ef f = h O0 q eit r + . . ., this pole at r* = —kj (/ioo^eff) hits the horizon and disap- 
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FIG. 9: The WKB estimate of the number of bound states n in the AdS-RN Schrodinger potential 
for z+ with mL = 10. The WKB approximation only applies to large values of the charge q = 45 
(Red), q = 50 (Blue), q = 55 (Purple). Fig B. gives the associated values of the IR conformal 
dimension = A= m 2 + k 2 — %^ . Both figures are for the extremal AdS-RN background with 
/i = \/3,r+ = l,g F = 1. 



pears precisely when k = . In AdS-RN, however, where lim^i fi qelf = + (r — 1) + . . . , 
the pole at r^ N — 1 = ~^ q {k+ hits the horizon and disappears at k — — -^=. For negative 
values of k whose magnitude is less than \k\ < the potential is regular and bounded and 
can and does have zero-energy solutions. Fig. IV shows this disappearance of the pole for 
the AdS-RN potential. 

Counting solutions through WKB is also more complicated for AdS-RN. For 0(1) values 
of q there are only few Fermi surfaces and the WKB approximation does not apply. For 
large q it does, however. For completeness we show the results in Fig. [9j 



V. CONCLUSION AND DISCUSSION 



These electron star spectral function results answer two of the three questions raised in 
the introduction directly. 

• They show explicitly how the fermion wavefunctions in their own gravitating potential 
well are ordered despite the fact that they all have strictly vanishing energy: In a 
fermionic version of the UV-IR correspondence they are ordered inversely in k, with 
the "lowest" /first occupied state having the highest k and the "highest" /last occupied 
state having the lowest k. With the qualitative AdS/CFT understanding that scale 
corresponds to distance away from the interior, one can intuitively picture this as 
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literally filling geometrical shells of the electron star, with the outermost/highest/last 
shell at large radius corresponding to the wavefunction with lowest local chemical 
potential and hence lowest k. 

• The decrease of the number of bound states — the number of occupied wavefunctions in 
the electron star — as we decrease q eff = fi l l A \J~^r for a fixed electron star background 
extrapolates naturally to a limit where the number of bound states is unity. This 
extrapolation pushes the solution beyond its adiabatic regime of validity. In principle 
we know what the correct description in this limit is: it is the AdS Dirac Hair solution 
constructed in [3]. The dependence of the number of bound states on k/L therefore 
illustrates that the electron star and Dirac Hair solutions are two limiting cases of the 
gravitationally backreacted Fermi gas. 

With this knowledge we can schematically classify the groundstate solutions of AdS Einstein- 
Maxwell gravity minimally coupled to charged fermions at finite charge density. For large 
mass mL in units of the constituent charge q, the only solution is a charged AdS-Reissner- 
Nordstrom black hole. For a low enough mass-to-charge ratio, the black hole becomes 
unstable and develops hair. If in addition the total charge density Q is of the order of the 
microscopic charge q this hairy solution is the Dirac Hair configuration constructed in [I], 
whereas in the limit of large total charge density Q one can make an adiabatic Thomas 
Fermi approximation and arrives a la Tolman-Oppenheimer-Volkov at an electron star (Fig. 

To). 

Translating this solution space through the AdS / CFT correspondence one reads off that 
in the dual strongly coupled field theory, one remains in the critical state if the ratio of the 
scaling dimension to the charge A/q is too large. For a small enough value of this ratio, 
the critical state is unstable and forms a novel scaleful groundstate. The generic condensed 
matter expectation of a unique Fermi liquid is realized if the total charge density is of the 
same order as the constituent charge. Following [T2J H3] and [TJHES] the state for Q ^> q is 
some deconfined Fermi liquid. 

The gravity description of either limit has some deficiencies, most notably the lack of an 
electron star wavefunction at infinity and the unnatural restriction to Q = q for the Dirac 
Hair solution. A generic solution for Q > q with wavefunction tails extending to infinity as 
the Dirac hair would be a more precise holographic dual to the strongly interacting large 
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Quant. Critical 



FIG. 10: Schematic diagram of the different groundstate solutions of strongly coupled fermions 
implied by holography for fixed charge density Q. Here q is the constituent charge of the fermions 
and mL ~ A the mass/conformal scaling dimension of the fermionic operator. One has the 
gravitational electron star (ES)/Dirac Hair (DH) solution for large/small Q/q and small mL/q 
dual a deconfined Fermi liquid/regular Fermi liquid in the CFT. For mL/q ~ A/q large the 
groundstate remains the fermionic quantum critical state dual to AdS-RN. 

N Fermi system. Any CFT information can then be cleanly read off at the AdS boundary. 
A naive construction could be to superpose Dirac Hair onto the electron star; in principle 
one can achieve this solution by a next order Hartree-Fock or Local Density Approximation 
computation. 

This best-of-both-worlds generic solution ought to be the true holographic dual of the 
strongly interacting Fermi ground state. If one is able to answer convincingly how this 
system circumvents the wisdom that the groundstate of an interacting many-body system 
of fermions is a generic single quasiparticle Landau Fermi liquid, then one would truly 
have found a finite density Fermi system that does not refer at any stage to an underlying 
perturbative Fermi gas. 
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